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We investigate quantum interference effects with magnetic field (magnetofingerprints) in strongly 
disordered Au-nanobridges. The magnetofingerprints are unconventional because they are caused 
by the Zeeman effect, not by the Aharonov-Bohm effect. These spin-based magnetofingerprints are 
equivalent to the Ericson's fluctuations (the fluctuations in electron transmission probability with 
electron energy) . We present a model based on the Landauer-Buttiker formalism that describes the 
data. We show that the dephasing time t$(E, T) of electrons at temperature T and energy E above 
the Fermi level can be obtained from the correlation magnetic field. In samples with localization 
length comparable to sample size, h/r^(E, T) ~ E, for E ^> fcsT, which shows that the Fermi liquid 
description of electron transport breaks down at length scale comparable to the localization length. 
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Studies of electron transport through very disordered 
nanometer scale materials are important for development 
of nano-devices. For example, nanometer scale bridges 
between noble metals have recently been fabricated us- 
ing electrodeposition and electromigration of noble met- 
als. 0,Q)H)SIEIE0)I3 I n certain cases, these nano- 
bridges can be highly disordered. The disorder can cause 
Coulomb Blockade or strong suppression in the density 
of states at low temperatures. jj-Ellil In this paper, the 
disorder is defined to be strong if the sample conduc- 
tance at zero bias voltage and low temperatures is much 
smaller than the conductance at room temperature. 

In our recent papers, we have reported on Coulomb- 
blockade and universal conductance fluctuations in 
strongly-disordered Au nanobridges formed by electromi- 
gration at room temperature in high vacuum. 0, [jj, [^j 
The disorder depends on growth conditions such as par- 
tial water vapor pressure in high vacuum and the ap- 
plied voltage. The Drude mean free path in the nano- 
bridge material was estimated to be ~ O.OIAf, showing 
that the disorder was significant. The origin of the dis- 
order was attributed to intermixing between the metal 
and H 2 molecules. We hypothesized that the disor- 
der was granular, because Au did not alloy with H 2 (so 
the disorder could not be amorphous). 

Recently, granular disorder has been directly observed 
in Au-nanoelectrodes grown by electrochemistry. 0, E3 
In particular, an electrochemical growth process at cer- 
tain gate voltages created a granular material consist- 
ing of Au clusters with oxidized surface. Granular disor- 
der can cause characteristics that resemble single-electron 
transistors, and researchers were cautioned about inter- 
preting transistor effects in molecular-electronics experi- 
ments. ^3 

The goal of this paper is to explain quantum inter- 
ference effects in strongly disordered Au nano-bridges at 
low temperatures, using a model of granular disorder. 
This paper is a follow-up to our recent letter [3] in which 
we demonstrated phase coherent phenomena in strongly 
disordered Au nano-bridges. We reported quantum inter- 



ference effects at the Fermi level similar to universal con- 
ductance fluctuations in weakly disordered metals. 

In the first part of the paper, we describe experimental 
results. The main observation is a linear structure in 
conductance versus bias voltage and magnetic field at low 
temperature. The structure demonstrates that electron 
interference effects with magnetic field are spin-based, 
in contrast to the similar effects in weakly disordered 
metals, which are caused by the Aharonov-Bohm effect 
on the phase of the electron wavefunction. 0] 

Spin-based interference effects were explained by the 
granular disorder model. || In this paper we present an 
in-depth discussion of the granular nano-bridge model. 
We introduce diffusion and resistor network approaches 
to granular conduction, and demonstrate the equivalence 
of the approaches, analogous to the Einstein relation in 
homogeneous systems. 

We show that the criterion to observe spin based mag- 
netofingerprints is that the Drude mean free path (Id) 
be smaller than the Fermi wavelength (Xf)- In homo- 
geneous systems, Id < Xf would imply that the system 
is an insulator at zero temperature. 14] However, in a 
granular system, Id can be much smaller than Xf even 
in the metallic state, because I n is not the same as the 
electron scattering length. [Til llflj This is why we can 
observe electron transport over relatively long distances 
(> 25nm) at the Fermi level. We present numerical simu- 
lations of quantum interference effects in a regime where 
I d < Xf, using a diffusive description of hopping con- 
duction. 

Finally, we present bias- voltage dependence of conduc- 
tance fluctuations, which has not been discussed yet. We 
obtain bias voltage dependence of the correlation mag- 
netic field, and show that this dependence displays the 
energy dependence of the dephasing time. In samples 
with resistance ~ Rq — h/e 2 , we find that the energy 
uncertainty of electrons at energy E above the Fermi 
level is the same as E. This shows that the Fermi-Liquid 
description of electron transport breaks down when the 
sample size is comparable to the localization length. 
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FIG. 1: A. Sketch of sample fabrication. B. Au nano-bridge 
with resistance of lOfcfi. 

EXPERIMENT DESCRIPTION 

In this section we describe sample fabrication and 
demonstrate spin-based quantum interference effects 
with magnetic field at low temperatures. 

Sample fabrication 

Nano-bridge fabrication process was discussed in detail 
in our prior publications, ja, 0, Ho| In summary, samples 
were created by deposition of Au atoms in high vacuum 
(~ 10~ 7 Torr) over two bulk Au films separated by a 
~ 70nm slit, as sketched in Fig.^-A. A finite bias voltage 
is applied between the films and the resulting current is 
measured in situ. The deposition of Au is stopped as 
soon as a finite current is detected. In this way, only one 
contact is created. 

The shape and the resistance of the contacts depends 
on growth conditions, such as bias voltage and back- 
ground pressure. For example, we demonstrated that 
the resistance of the contact at high bias voltage can be 
increased reversibly by several orders of magnitude by in- 
creasing the partial water vapor pressure around the con- 
tact. In general, the growth process is quite compli- 
cated, because processes such as electromigration, surface 
atom diffusion, and intermixing with water molecules af- 
fect the contact shape and the disorder. 

The contacts selected in this work are grown at 10V 
using a limiting resistor of R = 20fcf2. The purpose of 
the limiting resistor is to reduce the current as soon as 
the contact is created, which prevents catastrophic elec- 
tromigration processes. 

One typical sample is shown in Fig.^-B. The resistance 
of the contact is approximately inversely proportional to 
contact apparent width w. Q Thus, the resistance of the 
contact material is approximately uniformly distributed. 
The disorder length scale (grain size) is much smaller 
than the contact size L ~ 25nm. 

In wide contacts (w > 50nw), R < Rq = h/e 2 and 
samples do not display Coulomb Blockade at our lowest 
temperature (lbmK). In narrow contacts (w < 50nm), 
R > Rq and samples exhibit Coulomb Blockade at low 



temperatures. |6| Thus, in samples where w ~ 50nm, the 
electron localization length is comparable to L. 

The resistivity of the material inside the nano-bridge 
was estimated to be p ~ 10 5 fiilcm, Q which is larger 
than the maximum metallic resistivity of 10 3 p£lcm. |l7| 
The corresponding Drude mean free path (Id), obtained 
from p = mvF /ne 2 lD, is Id ~ O.OlAf . 

Spin-Based Magnetofingerprints 

We present quantum interference effects with magnetic 
field in devices with resistance < Rq, because these de- 
vices do not exhibit Coulomb Blockade at low tempera- 
tures, which makes it possible to study conductance fluc- 
tuations at the Fermi level, analogous to conductance 
fluctuations in weakly disordered metals. We select 
two new devices with resistances slightly below Rq. 

The I-V curves display a conductance suppression at 
zero bias voltage and at low temperatures (zero-bias 
anomaly or ZBA). This suppression is a Coulomb Block- 
ade precursor. |Q$ Figs. 0-A and B display differential 
conductance G = dl/dV versus V in samples 1 and 2, 
respectively. G is measured by lock-in technique with an 
excitation voltage 3pV. The ZBA was explained in detail 
by the theory of broadened Coulomb Blockade in a single 
electron transistor. 0] 

Figs. 0-C and D display differential conductance ver- 
sus magnetic field of samples 1 and 2 at zero DC-bias 
at base temperature. The conductance exhibits repro- 
ducible fluctuations with field. The fluctuation ampli- 
tude close to the universal value (~ e 2 /h). The fluc- 
tuations are sensitive to thermal cycling to room tem- 
perature, which demonstrates sensitivity of the fluctua- 
tions to impurity configurations. Thus, the fluctuations 
are analogous to the well known magnetofingerprints in 
mesoscopic physics. 

The magnetofingerprints in our samples are caused by 
the Zeeman effect. This is demonstrated in Fig.|3 which 
represents G versus magnetic field and bias voltage. The 
figure displays a linear structure, which is highlighted by 
dashed lines of mathematical form eV ± 2psB = const. 
We investigated a set of 10 samples with resistances in 
range (5feO, 15fcO) at milli-Kelvin temperature. The lin- 
ear structure was found in every one of these samples. 
Since the slopes of the lines correspond to the Zccnian 
splitting, this confirms that the fluctuations with mag- 
netic field are spin-based. By contrast, conductance fluc- 
tuations in weakly disordered metals are based on the 
Aharonov-Bohm effect and G versus V and G versus B 
are parametrically inequivalent. 

In the next 2 sections of this paper, we will demon- 
strate that magnetofingerprints become spin based when 
the Drude mean free path becomes shorter than the 
Fermi wavelength. If the conducting material were ho- 
mogeneous, this condition would not be physical because 
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FIG. 2: A and B: Differential conductance versus bias voltage 
at T — 0.015A" in samples 1 and 2, respectively. C and D: 
Differential conductance versus magnetic field at T — O.OlSiC 
in samples 1 and 2, respectively. 



the Drude mean free path and the electron scattering 
length would be the same. But, if the system is not 
homogeneous, the Drude mean free path and the elec- 
tron scattering length will not be the same and Id can 
be <C Xf , even if the conducting material is a metal at 
T = 0. 



GRANULAR NANO-BRIDGE MODEL 

In this section we introduce a model of electron trans- 
port through the nano-bridge. The model assumes that 
the disorder is granular with the grain size much smaller 
than L. The grain diameter (D) is certainly smaller than 
our imaging resolution of roughly 5nm. 

As discussed in the introduction, amorphous disorder 
is ruled out because Au does not alloy with the impu- 
rities that are responsible for disorder (H2O). Granular 
disorder has been recently confirmed in Au samples cre- 
ated by electro-chemistry. E3 It has been shown that 
large electric fields near Au surfaces induce formation of 
Au nanoparticles of diameter of few nanometers, which 
assemble into a granular structure. 01 

Theoretically, granular metals are similar to Fermi- 
Liquids. They are also known as Granular Fermi Liq- 
uids. |20j It has long been recognized that the Drude 
mean free path of a granular metal Id can be much 
shorter than the electron Fermi wavelength, 

[3E3 be- 

cause Id and the electron scattering length I are not the 
same. 

Consider a granular bridge of length L ^> D in Fig.QJ- 
A. The resistance of the bridge can be obtained us- 
ing a 3D-resistor network model, because, by definition, 
the resistance between the grains (Rg) is larger than 
the resistance inside the grains (otherwise, the system 



0.2\ 




5 B(T) 



10 



FIG. 3: A and B: Differential conductance versus magnetic 
field and bias voltage at T — 0.015if in samples 1 and 2, 
respectively. 



would be homogeneous). The resistance of the network 
is R = RgDL/A, where A is the sample cross-section 
area. The resistivity is obtained from R = pL/A, so 
p = R G D. 

In homogeneous systems, one can apply the Drude for- 
mula for resistivity p = ^fj- , where I is the electron scat- 
tering length. In granular systems, however, the Drude 
formula is not valid. But, one can still apply the Drude 
formula p = ™i{ D , where l D will be referred to as the 
nominal Drude mean free path. In a resistor network, 



In = 



mvF 
ne 2 R G D' 



(1) 



Id is not the same as electron scattering length I (I ~ D 
for ballistic grains) . One can change / d without changing 
I, by changing Rq- 

We expect that the metal insulator transition occurs 
roughly when Rq ~ zRq, where z is the coordination 
number [z = 6 in our model). This condition is reason- 
able, because if the total resistance between a grain and 
its neighbors exceeds Rq , then electron transport will be 
blocked by Coulomb-blockade and the system will be an 
insulator. Theoretical value for Rq at the metal-insulator 
transition is similar to our estimate. |2l| 

A granular metal will satisfy the condition Id < Xf if 
the inter-grain resistance is in interval Rq^d~ < < 
zRq. At the metal-insulator transition, Id has a critical 
value Xp/(zD) 1 which can be several orders of magnitude 
shorter than Xf- 

In this paper, we investigate mesoscopic electron trans- 
port in this regime (Id < Xp)- We will show that Id plays 
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FIG. 4: A. Mesoscopic granular conductor connected between 
reservoirs S and R. B. Quasi- Fermi level versus position when 
S is a source. 



a crucial role in determining the origin of mesoscopic fluc- 
tuations with magnetic field. 

Next, we discuss electron diffusion through granular 
systems. We assume that the average electron dwell 
time on any given grain is to- After time t 3> t&, 
there are roughly ^ hops between neighboring grains. 
If we assume that the hops are uncorrelated, then elec- 
tron motion is diffusive, and the diffusion length is L — 

D</i^- = \ r^-t. Thus, the diffusion coefficient is 
D dl f = 

In a homogeneous metal, the diffusion coefficient is 
r)dif _ V pijs_ In a granular metal, this relationship 
is not valid. One can still apply this formula, D dl f = 
^F^di//3, but Idif will not be the same as I and it will be 
referred to as the nominal diffusion mean free path. It 
follows that 
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Next, we investigate the relationship between the nom- 
inal Drude mean free path and the nominal diffusion 
mean free path. Let us assume that Idif — Id- Then 
= D R - G .p e and after some algebra it can be shown 
that this is equivalent to 



td 



hRa 

sr q - 



(2) 



This is the well known formula for the dwell time of 
an electron on a quantum dot connected to leads via re- 
sistance Rq. Eq- El is undoubtedly correct, so it follows 
that Idif = Id- Thus, when examining electron trans- 
port over distances much larger than grain diameter, one 
can use a diffusion model with Id obtained from the 
resistor- network model(Eq. QJ. Formula [5] is essentially 
the Einstein relation between diffusion and conductivity 
of a granular metal. 



MAGNETOFINGERPRINTS WHEN l D < \ F 

In this section, we explain the spin-based origin of con- 
ductance fluctuations (G(B, V, T)) and discuss linear cor- 
relations in differential conductance with magnetic field 
B and bias voltage V. First we examine fluctuations in 
differential conductance at T = 0. Finite temperature 
effects arc discussed in the appendix. 

We will neglect the effects of spin-orbit scattering. In 
our prior publication, we demonstrated that spin-orbit 
scattering is suppressed because of the quantization of 
energy levels inside the grains. 



Magnetofingerprints as the Ericson's fluctuations. 

First we investigate conductance fluctuations with field 
at zero DC bias voltage. We find the conductance contri- 
bution from electrons with spin-up, G^(B, 0, 0). The dif- 
ferential conductance is given by the Landauer-Buttiker 
formula, 22] 



Gi(B,0,0) = -T(E,B), 



(3) 



where 5 is the level spacing inside the grains and Rq 
h/e 2 is the resistance quantum. 



where T(E, B) is the transmission probability of elec- 
trons with energy E in magnetic field B. By our con- 
vention, the energy in this equation is measured relative 
to the Fermi energy Ep, i. e. E — for electrons with 
energy Ep- 

The magnetic field dependence in T(E, B) originates 
from the Aharonov-Bohm effect on the phase of the elec- 
tron wavefunctions. We will neglect the Aharonov-Bohm 
effect, so that the dependence T(E, B) on the second ar- 
gument can be omitted, i.e. T{E,B) = T(E). This ap- 
proximation is valid if the applied magnetic flux through 
the sample is much smaller than the flux quantum. Later 
in the text, we discuss the validity of this assumption. 

At zero bias voltage, E is just the kinetic energy. The 
dependence of the transmission probability on E arises 
because electron wavelength decreases with kinetic en- 
ergy. This causes mesoscopic fluctuations known as the 
Ericson's fluctuations. The Ericson's fluctuations in 
mesoscopic conductors have not yet been measured di- 
rectly. 
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Without Zeeman splitting, E in Eq. [3] would be equal 
to 0. However, with Zeeman splitting included, the ki- 
netic energy of spin-up electrons at the Fermi level in- 
creases with magnetic field. In particular, at the Fermi 
level the sum of the kinetic energy and the Zeeman en- 
ergy equals Ep. So E = g^ B B/2, where g = 2 is the 
g-factor, and /is is the Bohr magnetron. Consequently, 
Eq. |21 becomes 

G t {B,0,0) = jT(fx B B). (4) 

Taking into account electrons with spin down, the linear 
(zero-bias) conductance becomes 

G(B,0,0) = j(T(ti B B)+T(-» B B)). (5) 

This equation demonstrates that magnetofingerprints 
are equivalent to the Ericson's fluctuations. This behav- 
ior is possible only because we neglected the Aharonov- 
Bohm effect. The correlation energy Ec is related to the 
correlation field B C - E c = HbB C - At T = and V = 0, 
Ec is given by the Thouless energy Exh — hD ^ > ■ |l3| 

Conductance Fluctuations with Bias Voltage and 
the Ericson's fluctuations 

Now we consider the case when bias voltage V is ap- 
plied between the reservoirs. We choose a symmetric 



where f(E) = 1/ '(exp(E /fc^T) + I) is the Fermi func- 
tion at T — and summation over a indicates sum over 
electron spin polarization. 

In the Landauer formalism, electron transport is hori- 
zontal, so the kinetic energy increases with position (x), 
because the quasi-Fermi level decreases with x. The 
quasi-fermi level is defined as usual, as the position of the 
Fermi level that the electrons would assume if they were 
in local equilibrium (this would be the case if electron- 
phonon relaxation rate were much larger than the elec- 
tron transport rate). 

In good metals, the quasi-fermi level is very close to 
the electrostatic potential energy e<j){x) 1 where <fii x ) is 
the electrostatic potential and e is the electron charge 
(e < 0). We assume that the electric field inside the 
sample is uniform, so e<p{x) = eV/2 — eVx/L for < x < 
L. 



bias, in which the electrostatic potential of reservoir S is 
V/2 and the electrostatic potential of reservoir L is —V/2. 
We assume that reservoir S is a source (V < 0). One 
would naively expect that the fluctuations of G(B,V, 0) 
with V are similar to the Ericson's fluctuations, because 
increasing V increases the kinetic energy of electrons that 
participate in transport. 



We will assume that the energy is conserved in elec- 
tron transport, i. e. the electron transport is horizontal, 
as sketched in Fig. 0J-B. At V ^ 0, the quasi-Fermi level 
inside the sample varies with position. As a result, the 
kinetic energy varies with x, so strictly speaking V fluctu- 
ations will not be equivalent to the Ericson's fluctuations. 



The transmission probability depends on both electron 
energy and bias voltage. In further discussion, T(E, V) is 
the probability of transmission from S to R for electrons 
with energy E at bias V. E is measured relative to the 
average Fermi level of the two reservoirs, as shown in 
Fig.HB. 



In absence of inelastic electron scattering, the current 
through the sample is given by the Landauer-Buttiker 
formula, |22^|: 



(6) 

I 

The transmission coefficient can be expressed as 

T(E,V) = J2 T m,n( E >V), (7) 

where the sum is taken over incoming channels (to) in 
the left reservoir and all the outgoing channels (n) in 
the right reservoir, and T m _ n (E,V) is the transmission 
probability from channel m into channel n. T m . n (E, V) 
can be expressed in terms of the probability amplitude 
t m ,n{E,V) for transfer from channel m to channel n, 
T m ,n{E, V) — \t m . n {E, V)\ 2 . Finally, the probability am- 
plitude is equal to the sum over probability amplitudes 
of various Feynman paths connecting the reservoirs, [13 

t m AE,V)~Y. e ^ F(E,V) > ( 8 ) 
p 

where P indicates an electron path and <t>p{E, V) is the 



I(B,V,0) 



£ / T(E + a f i B B,V)[f(E 

a=±l ' 



2 ' 



-f(E 



eV 



) dE. 
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phase accumulated by an electron injected from reservoir 
S at energy E along this path. 

In the semiclassical approximation, <pp(E,V) = 4>p + 
K~ x f F SE(t)dt, where SE(t) is the kinetic energy of an 
electron on path P at moment t, and <f>p = <fip(0, 0). So, 
5E(t) = E-e<j> [x{t)\ and 

i ,„ t r\ i<t>p-iEt P /h+ieV C P (l/2-x(t)/L)dt/h 

tm,n\E, V ) ~ e ' 

P 

This expression can be rewritten as 

t m , n (E, V)~Y1 e i( " p - iEtp ' n+ieVtp ^ p V n , (9) 
p 

where = [1/2 - x{t)/L] dt. 

The dependencies of T(E, V) on P and V are not 
equivalent, because phase contributions Etp/H and 



T,(P)eVtp/h are not directly proportional, since E(P) 
fluctuates among different paths. If P is chosen so that 
an electron spends roughly the same time at different lo- 
cations x, as sketched in Fig. 0JA, then E(P) s=s 0. If 
an electron spends significantly more time closer to one 
reservoir than to the other reservoir, which is sketched 
by path P' in Fig. HA, then E(P') « 1. 

We determined the statistical distribution of E(P) 
among different paths. To this end, we generated 
random-walks on 3D grids connecting two reservoirs. We 
find that the distribution of £ is a Gaussian with a proba- 
bility density /(£) = (2ira 2 )- 1/2 exp(-E 2 /(2cr 2 )), where 
a = 0.167. The value of a is important, because it will set 
the voltage range where V-fluctuations and the Ericson's 
fluctuations are equivalent. 

Next we obtain the differential conductance 
G(B, V, T = 0) = dI(B, V, T) /dV, 



G(B, V, 0) = £ |- fdE-T(E + an B B, V) (-f'(E - ^-) - f(E 



eV, 
2 ' 



m E _ f ) _ /(B + £„. (10 , 

Q — ±1 



Substituting the Fermi function at T = 0, we find 



e 2 ^ el/ lei ^ f eV/2 

G(B,V 1 T = Q) = — T(a^ B B + a'—,V) + ^-J2 / dE 

1 q,q'=±1 a=±r- el/ /2 



^dT(E + au B B, V) 
97 



(11) 



The Ericson's fluctuations are represented by the de- 
pendencies T(ael//2 + a' /ipB, V) on the first argument. 
If both B and V are varied simultaneously so that 
eV/2 ± n B B = const, then two out of the four terms 
in the sum on the right hand side of Ea.llllwill not fluc- 
tuate, which will lead to a linear structure in G(B, V,0) 
similar to our experimental observations. 

The characteristic scales for the dependence of 
T(aeV/2 + a'[i B B, V) on the first argument are 

\e\V c /2 = fi B Bc=E c = E Th . (12) 

The linear structure in G(B, V, 0) will be significant if 
we can neglect the dependence T(aeV/2 + a'u B B, V) on 
the second argument (V) and if we can neglect the depen- 
dence of the integral on the right hand side of Eq. ^2 on 
B and V. Now we show that these two dependencies can 
indeed be neglected, which explains the linear structure 
in our experiments. In particular, we show that the cor- 
relation voltage Vq for these two dependencies is much 
larger than Vc- 



The correlation energy (Eq) is roughly given by condi- 
tion Ectp/h = ir (Eq. El for a typical path P. Similarly, 
the characteristic voltage (Vq) for the second argument 
in T(E,V) is obtained from condition eV^atp/h = tt 
(Eq. El). It follows that V$ « 6P c /|e| = W C - The 
dependence T(an B B + a'eV/2, V) on the second argu- 
ment has three times larger correlation voltage than that 
in the first argument, explaining the pronounced linear 
structure in our experiments. 

In conclusion, the conductance fluctuations with bias 
voltage are equivalent to the Ericson's fluctuations and 
Ec = |e|Vc/2. The equivalence persists only in voltage 
range < 3Vc- 

The integral in Eq.^Jalso leads to fluctuations in con- 
ductance with B and V. It has been predicted by Larkin 
and and Khmelnitskii |3| that this integral causes an in- 
crease in amplitude of conductance fluctuations at bias 
voltages larger than Ec/e, by a factor of ~ \/eVc/Ec- 
The increase has been observed in Ag wires. |2a| 

Now we discuss the characteristic V and B of the in- 
tegral. The integration over E broadens the dependence 
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FIG. 5: A. Transmission probability T(E,V) B. Differential 
conductance G(B,V). 

of dT(E + afisB, V)/dV on magnetic field, thereby in- 
creasing the correlation field from Be to eVj \±b- Thus, 
the correlation voltage of the integral is set by the V- 
dependence in T(E,V), and this correlation voltage is 
Vq, six times larger than Vc- So, the integral will not 
obscure the linear structure. 

The dependencies of T(u\xbB + a'-^-, V) on the sec- 
ond argument and the integral are responsible for the 
asymmetry of conductance fluctuations around zero bias 
voltage. 

To investigate different terms in Eq. 1111 we make nu- 
merical simulations. We model our sample with a grid 
containing 100 x 100 x 100 grains, and generate 2000 ran- 
dom walks across the sample. We calculate the phases 
4>(E, V) along these random walks and obtain T(E, V) 
from the transmission amplitudes in Eq. |U| We assume 
that there is only one conducting channel, for simplicity, 
and sum the transmission amplitude over 2000 randomly 
selected paths P. T(E, V) is plotted in Fig. EJA. The 
pattern is streched along V direction, which shows that 
that correlation scale for E (Ec) is much smaller than 
correlation scale for eV (\e\Vg), in agreement with the 
analysis above. 

The differential conductance versus zeeman energy and 
bias voltage are plotted in Fig.^B. The linear structure 
is clearly significant and similar to our experimental ob- 
servations. Note that there is a finite voltage range over 
which the lines in B — V parameter space are resolved. 
This voltage range is roughly 3Vc- 



Absence of the Aharonov-Bohm Effect 

We show that if Id < Xf, then quantum interference 
between paths P and P 1 in Fig. 0J-A will be more sen- 
sitive to the Zeeman effect than to the Aharonov-Bohm 
effect, which explains the spin-origin of magnetofinger- 
prints in our samples. The characteristic magnetic field 
for the Aharonov-Bohm effect is given by the field for 
a flux quantum $o = h/2e over sample area, Bab = 



a 




1 0.1 - 



o.o I 

100 200 300|iV 

V 

FIG. 6: Rms conductance versus bias voltage in sample 1 at 
T = 0.015A" . 

®q/L 2 . Using Eq.|H we find B C /B AB ~ Id/Xf, so 
B c < B AB if l D < A F . 



DEPHASING EFFECTS 

In this section we discuss electron dephasing in our 
nano-bridges. Samples can be divided into two groups. 
In the first group, the conductance fluctuation amplitude 
and correlation scales Be and Vc have strong T and 
V dependence. For example, in sample 1 conductance 
fluctuation amplitude decreases rapidly with \V\. Fig. 
displays root-mean-square differential conductance ver- 
sus bias voltage at T = 0.015^ in sample 1. The linear 
structure is pronounced only near V = 0, which indicates 
the the correlation field increases with \V\. Eight out of 
ten samples belong to the first sample group. 

In the second group of samples, Be and Vc at V = 
are independent of T below an experimentally accessible 
temperature. In addition, base temperature values of 
Be and Vc are independent of |V| below a certain \V\. 
This indicates that thermal broadening and dephasing 
are weak at low T and \ V\, so correlation scales are set by 
the Thouless energy, Ec = \eVc/2\ = \ibBc — Etk = 
hD dl f f / L 2 . Sample 2 belongs to this group. Two out 
of ten samples belong to the second sample group. As 
discussed before, all ten samples have similar resistance. 
The origin of different sample groups will be discussed 
later in the section. 

In the remainder of this section, we focus on the first 
group of samples. The dephasing time (r^) can be ob- 
tained from the voltage dependence of the correlation 
field at base temperature. This method for measuring 

is different from that in standard mesoscopic physics 
of weakly disordered metals, where can be measured 
from the temperature dependence of Be- |26| 

In our nano-bridges, magnetofingerprints are spin- 
based and fJ.Bc is equal to the correlation energy. In 
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Appendix 1 we find that if thermal broadening and de- 
phasing are significant, then 



HbB c (V,T) a 2M B T + h/r4E,T), 



(13) 



where E = \eV/2\. (The factor of 2.9 is valid in a narrow 
temperature range E Th < 2.9k B T < 3E Th ). t${E,T) is 
the dephasing time of electrons at energy E above the 
Fermi level at temperature T. 

It is not possible to obtain the temperature depen- 
dence of the dephasing time from the temperature de- 
pendence Bc(0, T), because Bc(0, T) is affected by both 
thermal broadening and dephasing. But, if T is fixed, the 
energy dependence of the dephasing time can be obtained 
from Eq. [TTTJ 

To increase the statistics, we thermally cycled sam- 
ple 1 between room temperature and base temperature 
8 times. Fig. shows the linear structure in 4 successive 
cool-downs, indicating that the data obtained in differ- 
ent thermal cycles are statistically uncorrelated. Even 
though thermal cycling scrambles the linear structure, 
the average conductance does not change with thermal 
cycling. 

Bc(V,T) is calculated as the field at which the field 
autocorrelation function (C(B,V,T)) decreases by 50%, 
i.e. C(B C (V,T),V,T) = 0.5. The field autocorrelation 
function at temperature T is found as 



C , B T) = G(B + B ,V, T)G(B , V, T) - G(B , V, T) 

G 2 (B ,V,T)~G(B Q ,V,T) 2 

where over-line indicates averaging over magnetic field 
interval < B < 10T. Fig. |S}A displays the field auto- 
correlation function at V — at base temperature. 

Next, we average the correlation field among differ- 
ent thermal cycles. The resulting correlation field versus 
bias voltage at base temperature is displayed in Fig. [HJB. 
Clearly, Bc(V,T) increases with V. In voltage interval 
< V < 2Q0[iV, the voltage dependence is approxi- 
mately linear, and the best linear fit is 

B C (V,T) = 0.28 + 4724V, 

where the left hand side is in Tesla and V is in volt. 

As discussed above, the voltage dependence of the 
correlation field displays the energy dependence of the 
dephasing time. Substituting E — eV/2 for the aver- 
age electron energy injected from the Fermi level of the 
source, we find 



fi B B c (V, T) = S2^ieV + 0.98 • E. 



(14) 



At energy E ^> 32[ieV, the offset in this equation can 
be neglected. Similarly, at large energy [E 3> k B T), the 
offset in Eq. ^] can be neglected, so comparing Eq. ^] 
and Eq. EI we find 



rt(E,T) 



0.98E. 



(15) 



(2) 



(3) 



(4) 




B(T) 
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FIG. 7: (l)-(4): Differential conductance versus magnetic 
field and bias voltage in sample 1 at T = 0.015^, obtained in 
four different thermal cycles. 



It follows that the correlation energy of an electron 
with energy E above the Fermi level is equal to E. Before 
we discuss the physical meaning of this observation, we 
examine the correlation scale at zero bias voltage. 

At E = eV/2 = 0, the correlation energy is E c (0, T) = 
HbBc(0,T) = 32fieV. This value is reasonable if we 
assume that the base electron temperature is T = 57mK. 
In particular, at V = 0, the average electronic energy is 
given by the full- width-half-maximum of the derivative of 
the Fermi function (3.5fc s T). So, E c (0,T) - 2M B T + 
0.98 • 3.5k B T = 6Ak B T. Substituting T = 57mK, we 
obtain Ec — 32[ieV. T = 57mK is a reasonable value 
for our base electron temperature, because the wires were 
filtered at the base phonon temperature (15mK). 

We believe that the equation h/T r p(E) = E demon- 
strates a break-down of the Fermi-Liquid description of 
electron transport in our samples. The selected sam- 
ples have conductance comparable to the conductance 
quantum. Coulomb Blockade in these samples is barely 
suppressed and the localization length is comparable to 
sample size. We expect that the Fermi-liquid picture in 
this regime begins to fail, which causes the energy uncer- 
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FIG. 8: A. Magnetic field autocorrelation function in sample 
1 at T — 0.015if . B. Correlation magnetic field versus bias 
voltage at T = 0.015^. 

tainty of a quasiparticle with energy E to be comparable 
to E. So, the correlation energy is comparable to the 
excitation energy. 

Now we discuss the variation of the Thouless energy 
among devices with a similar resistance. The Thouless 
energy is directly related to the electron diffusion time 
across the device, which implies that the electron diffu- 
sion time varies among different samples. The diffusion 
time is shorter than the thermal time (/i/fceT) in two 
out of ten samples with resistance slightly below the re- 
sistance quantum. 

The conductance-fluctuation amplitude is comparable 
to e 2 /h in all ten samples. This implies that in the eight 
samples where the Thouless energy is < k B T, the Thou- 
less energy cannot be <C k B T. If we assume the op- 
posite, then the conductance-fluctuation amplitude will 
be suppressed by \J ETh/k B T <C 1, |2j| contrary to the 
data. So, the diffusion time in our samples varies roughly 
within an order of magnitude. 

The fluctuations in the diffusion time among samples 
with similar resistance can be explained by the fluctu- 
ations in sample dimensions w and L. The fabrication 
process is based on electromigration and we do not have 
precise control over these dimensions. We found through 

I 



G(B,V,T) = - 




at T > 0. It follows that 

G(B, V,T) = J dx(-f'(x))G(x/n B + B, V, 0). (16) 

In this approximation, the conductance fluctuations at 
T > can be obtained by convolving the conductance 



scanning electron microscopy that the conductance is 
only roughly proportional to the apparent width w. The 
conductance per unit width varies within a factor of 3 
among different samples. Since the diffusion time scales 
as L 2 and the selected ten samples have similar resis- 
tance, the diffusion time can vary by factor of ~ 9, even 
if we assume that the resistivity does not vary among 
samples. This can explain the observed variation in the 
Thouless energy among samples. 

The fact that the fluctuation in the diffusion time 
among samples is within an order of magnitude indicates 
that our sample parameter control is also within an or- 
der of magnitude. We think that this sample control is 
quite good, in the light that the nanojunctions were self- 
created through electromigration. 



APPENDIX: DERIVATION OF THE 
CORRELATION FIELD AT FINITE 
TEMPERATURE 

At finite temperature, Ea. llQI still applies, but with / 
replaced with the Fermi function at temperature T. The 
formula takes into account thermal broadening effects, 
but it neglects dephasing. 

First we account for the thermal broadening. In Ea. llll 
which is valid at T = 0, we found that the correlation 
voltage of the integral is three times larger than that for 
the first term on the right hand side of the equation. 
Thus, to find the temperature dependence Vc(T) in the 
limit of small T, we only need to find the correlation 
voltage of the first integral in Eq. ^| The effect of the 
second integral in Eq. E3 on Vc (T) is weak if Vc (T) < 
3Vc(0) = 6E Th . 

If we neglect the integrals in Eas.lllland I1UI the con- 
ductance becomes 

G(B,V,0) = ^ T(a» B B + a'^,V) 

a, a' — ±1 

at T = and 



T(E + af i B B,V)f'(E + a'^) 

I 

fluctuations at T = with the derivative of the Fermi 
function. We reiterate that we neglected the conductance 
contributions with larger correlation scales. 

It follows that the field correlation function at temper- 
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ature T is proportional to 

dxdyf'(x)f'(y)C((x - V )j m + B, 0). 



where C is the field correlation function at T = 0. We 
can make an approximation C((x — y)/n B + ^>0) ~ 
S((x — y)/ Hb + B), because the thermal of /' (3.5k B T) 
is much larger the width of C((x — y)/fJ>B + B,0) (the 
Thouless energy). So, the field correlation function at 
temperature T becomes proportional to the autocorrela- 
tion function of the derivative of the Fermi function. The 
correlation field \x B Bc is the half- width- half-maximum of 
the autocorrelation function of /' (which is 2.9k B T), so 

Ms -Be = E c = 2.9k B T. 

In the presence of dephasing, the correlation energy 
is increased by energy uncertainty h/r,p. J^. So, the 
correlation field becomes 

fi B Bc = Ec = 2.9k B T+ —. 

In good Fermi liquids, one can neglect the dephasing 
term. 



CONCLUSION 

We show experimentally and theoretically that quan- 
tum interference effects with magnetic field in highly dis- 
ordered Au nanobridges are governed by the Zeeman 
effect. The Aharonov-Bohm effect is suppressed when 
the nominal Drude mean free path is smaller than the 
Fermi wavelength. Using the Landauer-Buttiker formal- 
ism, we show that both magnetofingerprints and conduc- 
tance fluctuations with bias voltage are equivalent to the 
Ericson's fluctuations. This permits measurements of the 
correlation energy from the correlation field. If sample 
size is comparable to the localization length, then the 
correlation energy of electrons at energy E 3> k B T will 
be ~ E, which shows that the Fermi liquid description 
breaks down at the localization length-scale. 
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David and Lucile Packard Foundation grant 2000-13874 
and the NSF grant DMR-0 102960. 
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